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asymmetric exclusion process with backward and forward site-ordered sequen- 
tial updates. This model, which was originally introduced with the random 



sequential update |11], describes a two-way traffic flow with a dynamic impu- 
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I. INTRODUCTION 

The one dimensional asymmetric simple exclusion process (ASEP) has been the subject 
of rigorous and intensive studies in recent years [0,0. A variety of phenomena can be 
modeled by the ASEP and its generalizations (see p|J^ and references therein). The model 
has a natural interpretation as a description of traffic flow ( |^ and references therein) 
and constitutes a basis for more realistic ones [^j. In traffic flow theories, the formation of 
traffic jams is one of the fundamental problems. Apart from their spontaneous formation |^, 
they can be also produced by hindrances such as road works or slow cars. Although these 
hindrances act locally, they may induce global and macroscopic effects on a system. This 
kind of behavior is one of the characteristic properties in non-equilibrium systems and has 
been studied in the context of driven lattice gases pj-p!^ . 



In the ASEP models, two kinds of impurities are discussed in the literature. The first 
one is ^^ dynamic impurities" , i.e., defective particles which jump with a rate lower than 



others p|,|T0|,|TT],|T3[ . In the traffic terminology, such moving defects can be visualized as 
slow cars on a road, which in certain situations can induce a phase transition from the 
free flow to the congested flow. The other kind of impurities is ^^static impurities" such as 
imperfect links where the hopping rate is lower than in other links |7|J9|JT^.T^-p!B|. Static 
impurities can also produce shocks in a system 0,^. For either type of impurities, a limited 
amount of exact results is available and most of them are for models with the random 
sequential update B JTO| , p!l|JT^ . For the fully parallel update, which is most realistic in traffic 



flow problems, exact solution are very rare 0,|^,0,|1^ and most studies instead utihze 
approximation methods or numerical approaches [0,0,|n]. 



Recently a two-way traffic model is introduced [^ where cars move forward in one lane 
and trucks move backward in the other lane, and both cars and trucks reduce their speeds 
when they approach each other. Within the matrix product ansatz (MPA) formalism, a 
modified version of this model is solved exactly in the particular case when there is only 
one truck in the system. This truck behaves as a dynamic impurity and induces a phase 
transition between the free flow and the congested flow of the cars. Various characteristics 
of the phase transition are examined. 

In this paper we study an exactly solvable traffic model with two types of ordered se- 



quential updates, which is identical to the model studied in Ref. [|TI[ except for the updating 
schemes. More precisely, the updating schemes that we consider are the backward and 
forward site-ordered sequential updating schemes in which one updates links of the chain 
sequentially. Alternatively one may use the so-called particle-ordered sequential updating 
scheme in which one sequentially updates the positions of the particles |]TB|. In this paper 
we restrict ourself to the site-ordered sequential updates. 

Presently it is of prime interest to determine whether distinct updating schemes can 



produce different befiaviors. Tfie implementation of tlie type of update is an essential part 
of the definition of a model and some characteristics of the model may change dramatically. 
The aim of this work is to investigate consequences of changing the updating scheme of the 
model. Then combined with the results in the random sequential updating scheme [11|, we 
examine similarities and differences between the results in different types of update. Our 
approach utilizes a mapping between the quadratic algebras of the ordered and random 



sequential updates, which has been initiated in the context of the one-species ASEP ||2T|j22 |. 
This paper has the following organization. In section II we define the model, construct 
our MPA with both backward and forward sequential updates, and obtain the relevant 
quadratic algebras. Section III presents the expressions of the average velocities of cars 
and truck, their thermodynamic limits, and a comparison with the corresponding results in 



Ref. [^. In section IV, we consider density profiles of cars and compute the probability 
to find a car at a distance x from the truck, which in the high density phase appears as a 
shock. We discuss the ambiguity in specifying the density profile, which is related with the 
nature of the updates, and introduce a new definition of the averages to avoid the ambiguity. 
Section V is devoted to the density-density correlation function. In section VI, we study the 
model in the presence of two trucks and evaluate the probability of their distance being R . 
The paper ends with some concluding remarks in section VII. 

II. MODEL DEFINITIONS AND MATRIX PRODUCT ANSATZ 

A. Two-way traffic model w^ith random sequential update 

Here we describe briefly the two-way traffic model with random sequential update (RSU) 
introduced in |0. Consider two parallel one dimensional chains, each with A^ sites. The 
periodic boundary condition applies to each chain. There are M cars and K trucks in the 
first and the second chain, respectively. Cars move to the right and trucks move to the left. 
The state of the system is characterized by two sets of occupation numbers (ri, T2,- ■ ■ , t^) 
and (cTi, o"2, ■ ■ ■ , <7n) for the first and the second lanes. If the site i of the car lane is occupied 
by a car, Xj = 1 and zero otherwise. Similarly ctj = 1 if the site i of the truck lane is occupied 
by a truck and ctj = if the site is empty. In an infinitesimal time interval dt, a car (truck) 
hops to its right (left) empty site with the probability dt {'jdt) if there is no truck (car) in 
front, and with the probability reduced by a factor (3 otherwise. More explicitly one has 

1 if o"j+i = 

(r„ r,+i) = (l, 0)^(0,1) with rate {l , , , . , ^ (1) 

— it (Jj+i = 1 (truck m front) 



7 if Ti = 

(cTi, ai+i) = (0, 1)^(1,0) with rate {j (2) 

— if Tj = 1 (car in front) . 

The reduction factor (3, which ranges between one and infinity, is related to the width of 
roads: (3 = 1 corresponds to a very wide road or a highway with a lane divider and (3 = oo 
corresponds to a one lane road. Simulations with finite densities of cars and trucks show [^ 



that in the steady state, the average velocities of cars and trucks decrease smoothly with 
increasing f3. As an interesting limiting case, situations with a single truck is considered 
while the density of cars is kept finite. For this particular case, simulations suggest that for 
a given density n of cars, there exists a density-dependent critical value f3c, below which the 
average velocity of cars remains constant and above which the average velocity decreases 
linearly with respect to 1 — 4, a measure of the road narrowness. For (3 > (3c, the simulation 
also finds the phase segregation into high (traffic jam) and low (free flow) density regions. 

To investigate the characteristics of this single truck case analytically, the above two- 
lane model has been modifled to an exactly solvable one. If one forbids a car and a truck to 
occupy two parallel site i simultaneously, one can describe conflgurations with a single set 
of occupation numbers {xj} where r^ = (empty site), 1 (occupied by a car), or 2 (occupied 
by a truck). The following rules describe the modifled dynamics: 



(3) 



This model is equivalent to a two-species ASEP and can be solved exactly by the method 
of the matrix product state (MPS). The steady state weight Pg of a given conflguration 
(ti, T2, ■ ■ ■ ? ''"A') is proportional to the trace of the normal product of some matrices: 

Psin, r2, ■ ■ ■ , r^) ~ Tr(XiX2 ■ ■ ■ X^) (4) 

where 

' D for Ti = 1 

X,= I 5 

and these matrices satisfy the quadratic algebra 

DE = D + E, aAE = A, (3D A = A {a = (3-f) . (6) 

With the help of the MPS method, many important characteristics such as average velocities. 
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density proflles, and k point correlation functions have been obtained exactly in Ref. [11 



B. Two-way traffic model with ordered sequential updates 

In the ordered sequential updating (OSU) schemes, time is discrete and the update occurs 
in an ordered way. To describe the dynamics in this scheme, it is convenient to introduce a 
Hilbert space that is spanned by the orthonormal ket vectors |{t}) = |ri) \t2) ®---® \'Tn)- 
The state of a system at the j-th time step can be represented by a ket vector |P, j) 

l^,j) = E^(M;j)IW) (7) 

conf 

where P{{t};j) is the weight of the configuration {r} at the j'-th time step and X^conf 
denotes the summation over all possible configurations. The state at the next time step is 
then obtained from |P, j) by applying a transfer matrix 

|P,j + l)=T|P,j) (8) 

where the transfer matrix T takes a different form depending on the precise nature of the 
updates. In the backward sequential updating (BSU) scheme, for example, the transfer 
matrix T becomes T^ where 

TV- = T]\jiTi2 ■ ■ ■ Tn-2,n~iTn-i,n (9) 

and in the forward sequential updating (FSU) scheme, it becomes T_^ where 

T^ = Tn,n-iTn-i,n-2 ■ ■ ■ T2^iTi^N ■ (10) 

Each element in the products is defined by 



-i,i+l 



Ti+i,i = l(g)---®l®T®l(g)---®l . (11) 



j-l N-i-l 



Here 1 is the identity matrix acting on a local ket vector |r) and the local transfer matrix T 
acts on the tensor product state of two local ket vectors, \t) ® \t'). In Eqs. (p|,p!OD, the site 
A^ is chosen as a starting point of the update for definiteness. 

The local transfer matrix T varies depending on possible exchange processes allowed in 



a model. As a straightforward generalization of the two-way traffic model in Ref. |TT|], we 
allow the following processes to occur in each discrete time step: 

(1, 0) — > (0, 1) with probability i] 

(0,2)^(2,0) with probability r/7 (12) 

(1,2)^(2,1) with probability ^ 



where < 77, 777, ^ < 1 < /3. For this two-way traffic model, off-diagonal matrix elements of 
T have nonzero values only for the following elements 

((0|®(1|)T(|1)®|0)) =v 

((2|®(0|)r(|0)®|2)) =r/7 (13) 

((2|®(1|)T(|1)®|2)) =j 

and the diagonal elements can be specified from the condition of the probability conservation, 
J2t3,t4 ((''"3 1 ® ('^4!) '^ (ki) ® 1^2)) = I5 which ensures that the sum of entries in each column 
is equal to one. 

C. Matrix product state 



Recently Krebs and Sandow |2^ showed that the steady state of any stochastic process 
with arbitrary nearest neighbor interactions, defined on an open one dimensional systems 
with RSU, can be always written as a matrix product state. This theorem was soon general- 



ized by Rajewsky and Schreckenberg p4| to stochastic processes with ordered sequential and 
sub-lattice parallel updates. For closed systems with the periodic boundary condition, on 
the other hand, such a theorem is not available at present and it is not yet clear how widely 
the MPA is applicable. Nevertheless, numerous recent studies H,p!0|, pJ] , ^ ,^5|-p7| have proven 
that the MPA can be a powerful investigation tool even for closed systems as well. Motivated 
by recent successes, we here assume that the stationary state of the two-way traffic model 
can be written in terms of the MPS: 

Psin, r2, ■ ■ ■ , tn) ~ Tr(XiX2 ■■■Xn) (14) 

where X, = D, E, or A {Xn = D, E, or A) depending on Xj (tn)- Note that the matrices at 
the site A^ are different from those at other sitesQ, which stems from the special role of the 
site A^ as a starting point of the update. The necessity of the hatted matrices in the OSU 
scheme is clearly explained in Ref. p2| . 

By introducing two ket vectors \U) and \U), 

\U) = A\0) + D\l) + E\2) 

\U) = A\0) + D\l) + E\2) , (15) 



^ The hatted-matrix in Eq. (^) breaks the translational invariance of the problem. As a result, 
the steady state weights are not completely fixed by the relative distances between cars and trucks. 
This situation is in contrast to Ref. [^I| where the weights depend on the relative distances only 
due to the translational invariance [see (^)]. 



where the coefficients A,D,E,A,D, and E are matrices defined on a different auxihary 
space, we can write ( |T^ ) formally as 

1^.) = 2^N^ M) Tr (|f/) ® If/) ® ■ ■ • ® \U)) (16) 

in which the normalization constant Z{N, M) is the sum of all weights 

Z(iV,M) = ^Tr(Xi---X^) (17) 

conf 

and the trace is taken only over the normal product of the matrices A, D, E, A, D, E. 
By definition, \Ps) should be stationary under the action of the transfer matrix T 

T\Ps) = \Ps) . 
In the BSU scheme, the stationarity is guaranteed if the following relation holds: 

Ti\U)(^\U)) = \U)®\U) (18) 

which simply implies that upon the action of T^, the ^^defecf \U) is transferred backward 
through the chain and returns back to the site A^, its initial position. Eq. (|l^) then leads to 
the following quadratic algebra: 

[A,A] = 0, [D,D] = 0, [E,E] = 
AD + r]DA = AD, {1-t])dA = DA 

r]-fAE + EA = EA, (1 - r]-f)AE = AE (19) 

^DE + ED = ED, {1-^)DE = DE 

In the FSU scheme, on the other hand, the stationarity requires 

T{\U) ® \U)) = \U) ® \U) (20) 

which implies that the ^^defecf \IJ) is now transferred in the forward direction. Eq. (|20|) 
leads to a quadratic algebra, which is identical to (p^Of ) upon the replacements A ■^^ A,D ^^ 

D,E ^ E. 

D. Mapping onto the RSU algebra 

In order to calculate physical quantities using the MPS, it is in principle necessary to find 
a representation of the matrices that satisfies the relevant quadratic algebras for the BSU 
and FSU schemes. In some cases, however, this job can be avoided and the quadratic algebra 

6 



itself is sufficient for calculations ||,|TT|,0,0,0,^ . Also recent studies on ASEP [0,^ 



and multi-species ASEP p8H30| have demonstrated that the algebra for a stochastic process 
in an OSU scheme can be mapped onto the algebra for a related stochastic process in the 
RSU scheme, for which an explicit representation of the matrices are known. Here we show 
that this mapping holds for the two-way traffic model as well. We ffist assume 

A = A + a, D = D + d, E = E + e (21) 

where a, d and e are real numbers. In the BSU scheme, it can be verified that the following 
choice 

(22) 



(23) 



Note that the algebra (23) is identical to the RSU algebra (P) except for the renormalization 
of a to a. 

In the FSU scheme, one may use the mapping to the BSU algebra by simply interchanging 
the roles of A,D,E with A,D,E [see Eqs. (p^SJpOl)]. However it turns out that it is more 



a = 0, d=-j3, 


e- ^ 




maps the algebra (pjj) to 






DE = D + E, aAE = 


= A, PDA = 


--A 


where 






_ a(/3 - v) 


R- R 





convenient for the subsequent analysis to use a direct mapping to the RSU scheme. We 
again assume Eq. (|2T|) and choose 

a = 0, d=-^, e = -|. (25) 

Then the FSU algebra is mapped to Eq. (p3D with 



a = a, /3 = ^-^ . (26) 

1 — ?7 

Note that now (3 is renormalized to (3. 

III. AVERAGE VELOCITIES 

In this section we consider the special case where there is only a single truck in the system 
and evaluate the average velocities of cars and the truck. 



A. BSU scheme 

The movement of cars from a site, k for example, to its neighboring site A; + 1 is achieved 
by Tfc^fc_|_i. Then the average current of cars {J^"^) between the two sites can be written as 
follows: 

-, ( N,M N,M 1 

«")=Z(]^ "S,^< ^ D^--)^IJ:M ^ DE-)\ (27) 

^ ' ^ I conf fc_i sites ^ conf fc_i sjtes ) 

where the specifications N and M in the summations represent that the matrix products 
(■ ■ ■ DA ■ ■ ■) and (■ ■ ■ DE ■ ■ ■) contain A^ matrices and M of them denote cars. Below, these 
specifications will be omitted in usual situations, that is, when they are N and M. In 
Eq. (0), it has been taken into account that before the action of T^^k+i, the steady state is 
acted by Tk+i^k+2 ■ ■ -Tn-i^n, which shifts the hat from the site A^ to A; + 1. More rigorous 
derivation is given in Appendix A. In a similar way, the average current of a truck {J^'^'^^) 
between the site k and k + 1 reads 

( Jf "^M _ _^ ^ J 777 y Tr( ■■■ A^ ■■■) + - y Tr( ■■■ D^ ■■■) i . (28) 

^ '^ ' Z(N, M)\ '^ ^, ^^^ P f ^^ I 

V ' ' \^ conf k-l ^ conf ^-1 ) 



It is instructive to briefly sketch the calculation of Eq. (|2^). Using the cyclic invariance 
and the fact that in the first term of Eq. (^), the single truck can locate anywhere on the 
remaining N — 2 sites, we find 

Note that the site index k has disappeared from the expression and thus the current is 
independent of k, as it should be from the conservation of the number of cars. Recalling 
that A = A and jSDA = A, the first term becomes 

N-2 N-l.M-l 
P 4=1 conf {_i 

In this double summation, each configuration {■ ■ ■ E) is counted [(A^ — 1) — (M — 1) — 1] 
times and thus it can be simplified further to 

1{N - M - l)Y{N - 1, M - 1) 



where 



P,Q 
Y{P,Q) = J2Tr{---E). (30) 

conf 



8 



Also recalling that E = E + e, the second term in Eq. (|29|) can be split into two pieces 

Ar-l,Af 
7] 7] 



lYniN,M) + e^ Y^ Tr(;^) (31) 

truck 



/5 ^ ' ^ P conf 



where 



P,Q 
YDiP,Q) = Y.M---DE). (32) 

conf 



Using (3DA = A and the cyclic invariance of the trace, the second piece in (RTl) can be 
simplified to 

where the binomial coefficient C^zi comes from counting the number of configurations that 
satisfy the given specifications. Later it turns out that the factor Tr(y4^~^^'^) cancels out 
from all expressions for physical quantities including (J^^^) and {Jl^^'^^)- Below we thus 
set Tr(A^~^^~^) = 1 for convenience. We also mention that due to the cancellation, all 
calculations in this paper can be performed without resorting to an explicit representation 
of the matrices. 

To evaluate the denominator of Eq. (P7|), we proceed with the definition ( [T7| ) 

Z{N, M) = Y: Tr(- . . E) + ^ Tr(- ■ • ^) + ^ Tr(- ■ ■ A) . (33) 

conf conf conf 

In a similar evaluation procedure as above, the first term on the RHS yields 

J:Tt{...E) = Y{N,M) + J^C^_,. (34) 

conf P 

Also the sum of the second and the third terms on the RHS becomes 

I:T.\m^D...E) + Tt{^A---E)\ , (35) 

4=1 conf [ j_i j_i J 

which results in 

{N - 1)Y{N, M) + d{N - l)Y{N - 1, M - 1) . (36) 

Combining all calculations given above, one obtains the exact expression for {J^^^) for 
arbitrary A^ and M: 

(4-)-^(^{|(iV-M-l)r(iV-l,M-l) 

+ |FMiV,M)+e|^C--} . (37) 

9 



For a complete evaluation of Eq. (^7]), one needs Y{N,M) and Yo{N,M). Exactly same 
quantities are calculated in the RSU scheme [|lT|. Taking into account the renormalization 
of a and f3, one immediately obtains 



Yd{N,M) 



^,RM 



a(3 



a 



1-/3 






a + (3-1 
/3-1 



/(iV, M) 



{31 



Y{N,M) = Yo{N,M) + ^^C'^_, , 



1 



M 



j(iv, M) = 5: ^p^c^'ii 



k ■ 



(39) 



fc=i 



To study the implications of Eq. (|3^), we consider the thermodynamic limit N, M ^>- oo 
while n = M/N is kept fixed. As demonstrated in Ref. [0, the thermodynamic limit is 
governed by I{N,M). Using the method of steepest descent, it can be verified that for 
n/3> 1 

37V-1 

I(N, M) ^ ^^— , 

^ ' {(3- l)N~M-i 

which is exponentially larger than all other terms and hence dominates the thermodynamic 
behavior of all quantities. For njS < 1, on the other hand, 

/(iV,M)^!^fcl!c«, 

1 — nfj 

which is of the same order as other terms. A straightforward calculation then shows that 
(^^car) = ^(-/rO is given by 

V 
— {I — II j 11 np ^ i 

m 

(40) 



(^^car) 



1 — rjn 
7] 1 — n 



1 - n) if n/3 < 1 
if n(3>l 



(3 — 7] n 
The average current of the single truck [Eq. (pSl)] can be evaluated in a similar manner: 



(4™^') 



1 



Z(iV, M) 

7] 



VI 



Y{N,M)-Yn{N,M) + J-C'J_, 



P 



M 



+ 



P 



Yo{N,M) + ^C^I,' 



e 



(41) 



In the thermodynamic limit, the average velocity (ftruck) = N{J^ ) becomes 



( 'Struck/ 



Tj {ail — 7i) + ea[7i + a(l — n)]} (1 — 7i(3) + n{a + (3 — n(3) . 



P 



{l-7]n) {l-7i){l-7i(3)+7i{a + P-Ti(3) 



V 



IP-V 



if np<l 
if 71(3 >1 



(42) 
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Figures 1 and 2 show the behaviors of (wcar) and (f truck) as a function of 1 — 4, the 
narrowness of the road, while the values of ^,7(= f) and n are fixed. The behavior of 
(■^car) changes clearly at the transition point 1 — -k- = 1 — n. Below the transition point, 
(■^car) is constant, implying that the interaction with the truck causes negligible changes 
on the flow of cars. Above the transition point, on the other hand, (fear) begins to drop 
suddenly generating a cusp at the transition point. Thus a local dynamic ^Hmpurity " , truck, 
results in global effects. Here the decrease of (fear) follows a hyperbolic curve, (ftruck) also 
changes its behavior at the transition point. Below the transition point, one finds a rather 
smooth decrease in (ftruck) and above the transition point, one again observes the hyperbolic 
decrease. 

B. FSU scheme 

In the FSU scheme, the average current of cars (J^^"^) between the site k and A; + 1 reads 

(Jr} = J^^{vJ:M^DA■■■)^'lY.T'(^DE■■■)\ . (43) 

V ' ^ \^ conf fc_i -^ conf j^_i ) 

In comparison with Eq. ( pT] ) for the BSU scheme, it should be noticed that the hats now 
appear at the site k instead of A; + 1 since the action of T^+i ^ occurs after the steady state is 
acted by Tk^k-i ■ ■ ■T2,iTi^n- This expression can be derived in a more rigorous way following 
a similar procedure as sketched in Appendix A. Using the quadratic algebra of the matrices, 
this expression can be reduced to 

<•''"> = Z(7P7) {" (i + <<) (iv - M - i)r (iv - 1. Af - 1) 

+1 [Yr,{N, M) + dY(N - 1, M - 1))| . (44) 

Note that the k dependence has disappeared. Here it should be understood that a,P,d,e 
now have the values given in Eqs. (|2^j2^). The same understanding is also needed for 
Y{N, M), Y£){N, M), Z{N, M) whose explicit expressions are given in the preceding subsec- 
tion. 

Similar shifts of the hats occur for the average current of the truck as well and one finds 

w> = z(i^{'"5;^(^-^^--)+i2:^(^^^--)} • <^^) 

V ' > \^ conf k_i I-' conf fc_i ) 

which is equal to 



11 



^^truck^ 



Z{N, M) 



7]^[Y{N,M)-Yd{N,M)] 



+ ^ [Yn{N, M) + dY{N - 1, M - 1)] 



(46) 



Note again that the k dependence has disappeared. 

Using the relations (fear) = ^{Jk^^) and (ftruck) = N{Jjf'^^^), and taking the thermody- 
namic hmit, one reaches the following results 



(Vc^t) 





r]{l-n) 


1 

V 


-ri{l-n) 

(l-n) 



. (3 n 



if n(3<l 
if n/3 > 1 



(47) 



("Struck) 



7] {a + dn(3){l - n){l - n(3) + {I + dn(3)n{a + (3 - n(3) 

T 

1 
1/3 



l + dnl3) {I - n){l - nl3) + n{a + [3 - nl3) 



if n(3<l 
if n/3 > 1 . 



(4J 



Figures 3 and 4 show (fear) and (ftruck) as a function of the narrowness 1 — ^ while 
?7,7(= ^), and n are fixed. One again finds that the average velocities have cusps at the 
transition point. Above the critical narrowness, both (fear) and (ftruck) decrease linearly 
with respect to the narrowness, which is in contrast to the hyperbolic decreases in the BSU 
scheme. Below the critical narrowness, (fear) is constant, similar to the result in the BSU 
scheme. But its value is different from the corresponding one in the BSU scheme. 

C. Comparison with the RSU results 



In the thermodynamic limit, (fear)Rsu and (ftruck)Rsu become HIT 

('i^car)RSU = 



1 — n if n(3 < 1 

lil-n) 

it np > 1 



(3 n 



(49) 



("Struck) 



RSU 



' 1 a(l - n){l - n(3) + n{a + (3- n(3) 
^ (1 - n)(l - nP) +n{a + P- np) 
1 



if np<l 
if nP>l . 



(50) 



We now compare the results. In all updating schemes considered, a phase transition 
occurs at a critical value of the narrowness. In the RSU and BSU schemes, the critical value 
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is 1 — n while in the FSU scheme, it is \Z^l[Sn) ■ Note that the transition point can be 
significantly lower in the FSU scheme when rj ^ 1. Below the critical narrowness, (fear) 
is independent of (3 in all three schemes (with diff'erent values in each scheme) and above 
it, (fear) decreases linearly in the RSU and FSU schemes but hyperbolically in the BSU 
scheme, (f truck)? on the other hand, is not constant even below the critical narrowness in all 
three schemes and varies quite smoothly with respect to the narrowness. Above the critical 
narrowness, (ftruck) decreases linearly in the RSU and FSU schemes but hyperbolically in 
the BSU scheme. 

In order to illustrate the origin of these differences, it is instructive to consider the 
behavior of (ftruck) in the vanishing car density limit n ^ 0. In this limit, (ftruck) in the 
RSU and FSU schemes approaches 7 and 777, respectively while it approaches j^ in the 
BSU scheme. This difference can be explained in the following way. In the absence of any 
car, the single truck in the RSU scheme either hops with the probability 'ydt or stays at the 
present site with the probability {l—'jdt) during a time interval dt. Hence its velocity is equal 
to its hopping rate 7. The situation is similar in the FSU scheme. Since the truck moves in 
the opposite direction of the update, it either hops one site ahead with the probabihty 777 
or stays with the probability (1 — 777). Therefore, the average velocity reads 

, , average distance x (1 — 777) + 1 x (777) 

(^^truck) = — r — rr- r~ = ^ = ^^ • 

number of time step 1 

The situation changes drastically in the BSU scheme. In this case, the direction of the 
update and the direction of the truck movement are identical and accordingly the truck can 
be transferred by large distances in a single time step. Recalling that each hopping occurs 
with the probabihties 7/7 and the probability of stopping is (1 — 777), (ftruck) can be cast into 
the following form 

X (1 - 777) + 777(1 - 777) + 2(777)^(1 - r/7) + 3(7/7)^(1 - 777) + ■ ■ ■ 



which is simplified to 7 

We next examine the fundamental diagrams (the relation between the car current and 
the car density) in three updating schemes. The average current of cars is simply related 
to the average velocity via (Jcar) = n{vca.r)- Using the thermodynamic behaviors of (fear) 
[Eqs. (0,^,^)], one can determine the thermodynamic behaviors of the current as a func- 
tion of the density. Figure 5 shows these behaviors for some constant values of f3 and 7/ in 
different updating schemes. In all three types of update, the single truck does not affect 
the current-density diagram before the transition point whereas it affects the system after 
the transition point in a nontrivial manner: linear decrease of the current with increasing 
density. 
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IV. DENSITY PROFILE 



A. Density profile in the RSU scheme 



Here we first summarize the results in Ref. 111. In the RSU scheme, it can be assumed 



without loss of generality that the single truck is at a particular site, for example at the site 
N of the chain due to the cyclic invariance of the problem. Thus the probability to find a 
car at the distance x from the truck can be written via MPS as follows 



{n{x)) 



RSU 



EconiM-'-D-^E) 



Sconf Tr( 



E) 



Using the matrix algebra (^), one finds 



(n(x)) 



RSU 



YRSv{N,M)aP 



M 



^M-1 
'-'N-3 



« r^M-2 ,a + l3-l 

'-'N-3 "I" 



P-1 



(3-1 



(51) 



X 



/Rsu(iV - 1, M - 1) + I3^-\I3 - l)/Rsu(iV - X, M - x)e{M > x) 



(52) 



where Y^sv{N,M) and Iksv{N,M) can be obtained from Y{N,M) and I{N,M) 
[Eqs. (p^,^)] by replacing a,/5 with a,/3, respectively and the factor 9{y > a;) is 1 if y > a; 
and otherwise. In the high density phase nf3 > 1, the single truck affects the system 
globally. The density profile is 



(n(x)) 



1 for f < /rsu 



n/3-l 
/3-1 



RSU 



/3 



otherwise . 



(53) 



Note that the system consists of two regions, a traffic jam region in front of the truck and a 
free flow region behind it. In the low density phase n(3 < 1, on the other hand, the presence 
of the truck has only local effects. In the thermodynamic limit, the car density becomes 



^^ a; RSU = n\l + f^ {npY 

1 — n + an 



(54) 



which shows that the disturbance by the truck decays exponentially with a characteristic 
length scale 



e 



RSU 



ln(n/3)|-^ . 



(55) 
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B. Density profile in the ordered sequential updating schemes 

As stated in Section 2, the ordered sequential updates break the cychc invariance. In 
what follows we show that due to the absence of the cyclic invariance, the probability to find 
a car at x sites in front of the truck varies depending on the truck location. Two updating 
schemes, BSU and FSU, will be considered simultaneously since same expressions apply to 
both schemes. 

We first consider the case where the truck is at the site N. The conditional probability 
to find a car at x sites in front of the truck reads 

X 

ProbfAT -l-x = carliV = truck) = ^confTr(---D ■■■E) ^ 

EconfTr(---E) 

Here it is helpful to define a quantity K{x, N, M) by 

/f(..iv.M) = g-'^-'-;-^^^f 

Ec„„tTr( E) 

which is exactly the same as Eq. (pi]). All properties of K{x,N,M) can be thus obtained 
from Eq. (^2]) by replacing a and (3 with a and (3, respectively. In terms of K{x, N, M), the 
conditional 

probability ( p6D becomes 



YiN,M)Kix,N,M) + ^C^ 



P^ 



M-l 

■'N-2 



Y{N,M) + ^C_ 



(3^ 



M 



(57) 



The thermodynamic limit can be investigated in a simple way. In the high density phase 
njS > 1, the second terms both in the numerator and the denominator in Eq. (|57|) are 
negligible compared to the first terms and Eq. (^Tf ) reduces to K{x, N, M). Thus one finds 



1 for ^ < / = "/^-^ 

Prob(A^ -x-1 = caT\N = truck) = <^ 1 (58) 

— otherwise 

P 

which is essentially identical to the result (p3|) in the RSU scheme except for the replacement 
of (3 by (3. In the low density phase n/S < 1, on the other hand, the second terms are 
comparable with the first terms and one obtains 

Prob(Ar - a; - 1 = car|Ar = truck) = n\l + (a + P - l)il - n) ^ ^^~^.J ^^^^ 

[ (1 — ra + an) + e(l — nf3)a J 

which shows the exponential decay of the disturbance with a length scale ^ = In (n(3] 
The coefficient of the exponential decaying term is different from the result in the RSU 
scheme. 
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Secondly we consider the case where the truck is at the site N — k (1 <k<N — x — 2). 
The conditional probability becomes 

X fc— 1 

ProbfiV -x-k-l = cailN -k = truck) = ^conf Tr(- ■ ■ Z^ ■■■ E ■ ■ ■ Xjv) ^ 

EconfTr(---E^X;v) 

fc-i 

The evaluation of Eq. ( |60| ) using the matrix algebra leads to 

F(Ar, M)i^(x, AT, M) + rfF(iV - 1, M - l)ir(x, iV - 1, M - 1) 

Y{N, M) + rfF(iV - 1, M - 1) ■ *^^"^'' 

Note that the conditional probability is independent of k. 

In the thermodynamic limit, (|6T|) can be greatly simplified to K{x, N, M) since 
K{x,N,M) and K{x,N — 1,M — 1) in the numerator become identical in this limit [see 
Eqs. ( |53| , |5^ ]. Thus the conditional probability Prob(A^ — x — k — 1 = car|A^ — k = truck) is 
the same as the corresponding results ( p3| , p4D in the high and the low density phases except 
for the renormalization of a and (3 to a and (3. 

In the third case, the truck is located at the site x + 1 and the car at the site A^. The 
conditional probability becomes 

X 

Prob(iV = carjx + 1 = truck) = ^conf Tr( • • • E- ■ -D) ^ ^g^) 

Econf Tr(^^-E ■ ■ ■ Xjv) 

X 

which is found to be 

Y{N, M)K{x, N, M) + dY{N - 1, M - 1) 



(63) 



Y{N,M) + dY{N -l,M -I) 

Note that the second terms both in the numerator and the denominator are comparable to 
the first terms in the thermodynamic limit. As a result, the conditional probability becomes 
noticeably different from K{x, N, M). One finds in the high density phase 

Prob(A^ = car|x + 1 = truck) = { I I + df3 (64) 

— otherwise 

I p l + d 

and in the low density phase 

X n(l + d(3) ( (a + 3-l)(l-n) . ~-sx] ^ ^ 

Prob A^ = carb + 1 = truck = ^ ^J \ 1 + ^ ^ (np) \ . 65 

^ 1 + dnp \ {l + dp)il-n + an)^ ^ j 

In contrast to the previous cases, the renormalization of a and (3 are not sufficient to account 
for deviations from the RSU results even in the high density phase. 
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Lastly we consider Prob(A^ — / = car|x — / + 1 = truck) (1 < / < a;), which takes the 
following form 

Y{N,M)K{x,N,M)+dY{N-l,M-l)K{x-l,N-l,M-l) 

Y{N, M) + dY{N - 1, M - 1) ' ^^^' 

Note that it is independent of /. The thermodynamic limit can be examined in a similar 
way. In the high density phase, K{x,N,M) and K{x — 1,N — 1,M — 1) are essentially 
identical and thus Eq. (^) reduces to Eq. (^) for the RSU scheme except for the trivial 
replacement of a and (3 with a and (3. In the low density phase, however, the difference 
between K{x, N, M) and K{x — 1, N — 1, M — 1) is not negligible and careful treatment is 
required to take care of the difference. This way, one obtains 

ProbfiV - / = carlx -1 + 1 = truck) = n |l + ^ + ^ J« + /^ - 1)(1 " ^) (^^v] (57) 

[ 1 + dn/3 l-n + an ^ ^ J 

We next discuss the origin of the absence of the cyclic invariance, which results in the 
four different cases. As stated above, the system loses the cyclic invariance due to the choice 
of a particular site, the site A^, as a starting point of the update. 

Thus by choosing the starting point in an even way, the cyclic invariance can be restored. 
One can, for example, define a cyclically invariant average ((•■■)) of an operator O in the 
following way0: 

1 ^ 
{iO)) = -j:{^\0\Ps,,) (68) 

^^ fc=i 

where {s\ = I]{r} ({''"} I ^"^^ \Ps,k)k is the stationary state of the transfer matrix 

rri rri rri rri rri rri 

-t «-,fc = -t k,k+l-'^k+l,k+2 ■ ■ ■ -t Af,!-' 1,2 ■ ■ ■ -t k-l,k 

in the BSU scheme or 

rri rri rri rri rri 

-^,k = -l-k,k-l'Lk-l,k~2 " " " J I.AtJ Af.AT-l " " " Jfc+l.fc 

in the FSU scheme with the site fc as a starting point of the update. Note that the site N 
now loses its special meaning and the cyclic invariance becomes evident in Eq. (|6^). It is 
also worth mentioning that the definition (^) is equivalent to taking an expectation value 
at each sub-step of an update, that is, after the action of Ti^i^i or Tj+i^/ instead of a single 
whole step of the update T^ or T_^, and taking the average over these expectation values. 



§This problem due to the absence of the cyclic invariance does not affect the calculation of the 
currents since the currents should be independent of the sites where their values are evaluated due 
to the conservation of the particle number. 
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Using the definition (|68|), we calculate the density ((n(x))) = ((5riv_i,2'^TAr_,_^_i,i)) 
/{{^TN-1,2)) (0 < / < A^ — 1), where the sites 0, — 1, — 2, ■ ■ ■ should be identified with the 
sites A^, A^ — 1, A^ — 2, ■ ■ -, respectively. The density {{n{x))) is, by definition, independent 
of the truck location A^ — / and thus one may choose / = 0. The denominator ((5t-jv,2)) is jj^ 
due to the cyclic invariance and the numerator (((5T-^,2'^r]v_j;_i,i)) becomes 



Af-x-3 



\ ' 'I com X /e=0 com i. x 



x-1 



conf X 1=0 conf i x—l—1 



(69) 



The cyclically invariant density then becomes 

1 



((n(x))) 



Z{N,M) 



NY{N, M)K{x, N, M) + -^C^~2 + dY{N - 1, M - 1) 



(3 



M 



+ d{N -X- 2)Y{N - 1, M - l)K{x, N -l,M -I) 
+ dxY{N - 1, M - l)ir(x - 1, A^ - 1, M - 1) 



(70) 



which, in the thermodynamic limit, reduces to the results identical to the thermodynamic 
behaviors of the second conditional probability Prob(A^ — x ~ k — 1 = car|A^ — k = truck) 
(1 < A; < A^-x-2). 

V. DENSITY-DENSITY CORRELATION FUNCTION 



Here we calculate the two-point equal time correlation function of the car density using 
the cyclically invariant average (p8|). Both BSU and FSU schemes are considered simulta- 
neously. In terms of the MPS, the density-density correlation function becomes (xi < X2 is 
assumed) 



{{n{xi)n{x2))) 



Z{N,M) 



X2 



^ Tr(- • • D ■ ■ ■ D ■ ■ • E) 



conf 



XI 



N-X2-3 



X2 



fc=0 conf xi k conf xi ]\[-x2-2 

+ '2 J2'^<---D^_:^E ^^ D^^Xm) + Y.M^^E ^^ D---D) 

k=0 conf xi N-X2-2 k conf xi N-X2-2 

xi-1 N-X2-2 

+ y: y^ri---E ^ ^ D.-- D,..x^] 

k=0 conf N-xi-2 k 
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After some algebra, it can be verified that 

{{nix,)nix2))) = {{nix2))) ~ fix,) (72) 

wliere f{xi) is given by 

f{x,) = ^^^^^^ i^j^C^li + dY{N~l,M-l)[l~K{x,,N-l,M-l)] 

+iNY{N - 1, M - 1) [1 - K{xi, N-1,M-1)] (73) 

+i{N -xi- 3)Y{N - 2, M - 2) [1 - K{xi, N - 2, M - 2)] 

+ixiY{N - 2, M - 2) [1 - K{xi - 1, A^ - 2, M - 2)] I . 
P ) 

Tlie connected part of the two-point correlation function, {{n{xi)n{x2)))c = 
{{n{xi)n{x2))) — ((n(xi)))((n(s2))), can be used to estimate the degree of correlation. In 
the thermodynamic limit, one finds 

{{n{x,)n{x2)))c = [{{n{x2))) - k][1 - {{n{x,)))] (74) 

where n = min(ra, i). In the high density phase, the connected part has non vanishing 



value only when both xi and X2 are within the region [Nl — yNA, Nl + V A^A] where 
A = y 2(3(1 — n)/((3 — 1). In the low density phase, the connected part has non vanishing 
value only when xi < X2 ~ ^. Thus one concludes that the correlation develops only in the 
region where the density variation occurs, which is identical to the conclusion in the RSU 
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VI. TWO TRUCKS AND BOUND STATE 

In this section we examine the system with two trucks and M cars. We determine the 
probability Q{R) that the distance between the trucks is R [0 < R < [N — l)/2]. In the 
RSU scheme, the cyclic invariance of the MPS allows one to set one of the trucks at the site 
A^ and one obtains 

f^RSu(i?)~ETr(---^^;^^) 

conf p> 

where the sum runs over all configurations with M cars and two trucks. Its thermodynamic 



limit is examined in Ref. |Tl|. In the low density phase n/3 < 1, it becomes 



r^ /T^\ . n(l ~ n)(a + (3 - l)(a - I) . ^,n ,„^, 

(1 — n + anj^ 
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which is maximal at i? = and decays exponentially with the same length scale ^rsu = 
I ln(?T,/5)|~^ as in the density profile. In the high density phase n(3 > 1, the probability 
decreases linearly with i? for < i? < A^rRsu, where trsu = niin(/Rsu, 1 - hsv) (^rsu < |), 
and remains constant for A^rRsu < -R < {N — l)/2. The relative ratios are: 

for rRsu = /rsu, — t^ ttt. — = 1 —5 (76) 

"RSu(O) «P 

and 

r , , ^RSviNr-nsv) , (a-l)(/3-l) r^sv (^^. 

for trsu = 1 - ^Rsu, — p; ttt. — = 1 7, ^ ■ 77 

i^RSu(O) a/3 1-rRsu 

One can interpret this as the formation of a weak bound state between the two trucks. 

Now we consider VL{R) in the BSU and FSU schemes. In terms of the matrix products, 
Vt{R) in both schemes can be expressed (up to a proper normalization constant) as 

Af-i?-3 

^ ^Tr(---^^^^^^X^) + ^Tr( ^^ ^^^^) 

fc=0 conf ji fc conf N-R-2 R 

R-1 

+ ^5]Tr(---E ^^ E^_^XN) + J2Tii^_^E ^^ E) . (78) 

A;=0 conf N-R-2 k conf R N-R-2 

It is useful to introduce a quantity W{N, M, R) which is defined by 

W{N, M,R) = Y,i:i{---E ■■■ E) . (79) 

conf ^ 

Then using the relation A = A,D = D + d,E = E + e and the cyclic invariance of the trace, 
the expression (|78|) can be written in terms of W{N, M, R) as follows: 

NW{N, M, R) + d{N -R- 2)W{N - 1, M - 1, i?) 

+dRW{N - 1, M - 1, i? - 1) + 2eY{N - 1, M) . (80) 

The thermodynamic limit can be investigated in a simple way using the fact that in 
the RSU scheme, W{N,M,R) is identical to QiisviR) up to a normalization factor. Then 
through the renormalization of a and (3, its R dependence in the BSU and FSU schemes can 
be obtained. Also the last term in the expression ( pOf ) is negligible compared to the first three 
terms. In the high density phase, one then finds that W{N, M, R), W{N — 1,M — 1,R), and 
W{N — 1,M — 1,R — 1) are all proportional to each other and thus Q{R) can be obtained 
from r2Rsu(-R) through the ^^renormalization" of a and (3. In the low density phase, the R 
dependence of W{N, M, R) appears only for i? ~ ^. It is then sufficient to examine the 
case i? ~ ^ <^ A^, where the first and the second terms are dominant and give the same 
R dependence. Hence Q{R) can be again obtained from Q^su{R) by replacing a and /? by 
their ^^renormalized^ values. 
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VII. CONCLUDING REMARKS 

We have investigated the characteristics of an exactly solvable two-way traffic model 
with the ordered sequential updates (OSU) and observed both qualitative and quantitative 
differences in the properties of model from the results obtained with the random sequential 
update |TT|. Our approach is based on the so-called matrix product formalism which allows 



analytic solutions. In the OSU schemes, the choice of a particular site as a starting point 
of the update breaks the translational invariance of the steady state measure, which is also 
evident in the form of the MPS. Thus an averaging over the different choices of the update 
starting point is necessary to restore the cyclic invariance to the system [see Eq. (^8]) and the 
following discussion]. Performing the cyclically invariant averaging, some characteristics in 
the thermodynamic limit, such as density proffie of cars, density-density correlation function, 
and truck-truck distance distribution Q{R), are obtained and it is found that the difference 
in the updating schemes can be taken into account simply by the proper renormalization 
of the parameters a and (3. However this is not the case with average velocities. Changing 
the update scheme affects velocities in a more complicated manner and the renormalization 
of the parameters is not sufficient to account for different behaviors of (fear) and (ftruck) 
in different updating schemes. Especially the dependence of (fear) and (ftruck) on the road 
narrowness 1 — 4 can vary qualitatively depending on the updating schemes. Behaviors of 
(fear) and (ftruck) in the FSU scheme (that is, when the update direction is parallel with the 
movement direction of the majority of vehicles) have more resemblance to the RSU results 
than those in the BSU scheme. In the FSU scheme, however, one observes a shift in the 
value of the critical narrowness from 1 — -n, to ^ i "" , which can be considerable if r? ^ 1. 
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APPENDIX A: DERIVATION OF THE EXPRESSION (27) 



Here we derive the expression (|27|) for the average current of cars in the BSU scheme. 
Similar procedures can be used to obtain average currents of the truck and also the average 
currents in the FSU scheme. The starting point is the continuity equation, which in discrete 
time dynamics takes the form 
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{nDj+i - {nT), = (Jr-i). - {Jn^ (Al) 

where (■ ■ ■)j represents the average at the time step j. In terms of the initial state of the 
system |-P, 0), the LHS becomes 

{s\nl^'T^Tl\P, 0) - {s\nl^'Ti\P, 0) (A2) 

where the bra vector {s\ is defined by 

(s| = E(^il®(^2|®---®(riv| • (A3) 

{r} 

The conservation of the probabihty ensures that {s\T^ = {s\ 0, which then allows one to 
rewrite ( [A2D as {s\[nl^^,T^]T^\P,0). Next we evaluate the commutator [nl^^,T^]. Using 
the relation 

nl""' = lOl® ••■(g)l®|l)(l| ® 1® 1® ••■® 1 , 



fc-l N-k 



one finds 



7^,1 ■ ■ ■ Tk-2,k-i ( vEk'-iE^' + -Ef^LiEf^' 1 Tk,k+i ■ ■ ■ Tn~i,n 
- Tn,i ■ ■ ■ Tk-i^k ( rjEk -£^fcVi + -n^l'^El'^i j Tk+i,k+2 ■ ■ ■ Tn-i,n (A4) 



where 



E]:^ = l(g)l(g)--- (g)l ®\i){j\ (g) 1 (g) 1 (g ) ■ ■ ■ (g) 1 . 

fc-l N-k 

Comparison with the RHS of Eq. (|A1|) (note that the two currents have different subscripts 
k — 1 and k) shows that each term in (|A4D should lead to the expression for {Jl^^j and 
{Jk^^)j, respectively. Then by taking the limit j —>■ oo and using the expression (p!6D for the 
steady state \Ps) = limj^oo T^\P, 0) , one finds 

'( ) 

xTi \U} «>\U} <S - ■ ■\U}<S>\U} <»\U} » - ■ ■ rS\U} (A5) 

v^ 1 ■ ) 

where the effects of T/^/+i on IP,) and (s| have been taken into account. Finally we use 

\t){j\{Ao\0) + A^\l) + A2\2)) = Aj\t) , 
which yields the expression (|27|). 
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FIGURES 

FIG. 1. Average velocity of cars in the BSU and RSU schemes for n = 0.3 and different values 
of 77 

FIG. 2. Average velocity of the truck in the BSU and RSU schemes for n = 0.3 and different 
values of t] 

FIG. 3. Average velocity of cars in the FSU and RSU schemes for n = 0.3 and different values 
of r] 

FIG. 4. Average velocity of the truck in the FSU and RSU schemes for n = 0.3 and different 
values of t] 

FIG. 5. Fundamental diagrams in three updating schemes for different values of (3 and rj (values 
are given in the figure) 
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